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I INTRODUCTION 

In  a series  of  observational  and  theoretical  studies  of 
the  air  flow  over  an  island  (Malkus  and  Bunker,  1902),  (Malkus  and 
Stern,  1953),  (Stern  and  Malkus,  1953) » the  connection  between  the 
turbulent  heating  in  the  lower  levels  of  the  atmosphere  and  the 
mean  velocity  perturbations  arising  from  this  heating  has  been 
brought  out.  These  mean  motions  may  appear  as  quite  sizeable 
perturbations  and  at  large  distances  from  the  energy  source,  as 
it  is  evidenced  by  the  regularly  spaced  cloud  streets  that  have 
been  observed  in  the  lee  of  small  flat  oceanic  Islands. 

However,  the  problem  is  of  more  general  Interest  than 
explaining  local  phenomena  due  to  heating.  It  is  hoped  that  if 
a satisfactory  model  of  the  small  scale  effects  can  be  obtained, 
then  tnis  may  eventually  lead  to  an  understanding  of  larger  scale 
and  more  important  meteorological  problems. 

A major  problem  In  the  previously  mentioned  theoretical 
studies  was  the  specification  of  the  gross  features  of  the  tur- 
bulent heating.  It  was  felt  that  the  conventional  eddy  conduction 
equation,  which  Is  frequently  used  to  get  a first  orientation  into 
many  problems  of  atmospheric  turbulence,  was  inadequate  In  the 
present  type  problem.  In  this  paper  we  shall  elaborate  on  this 
inadequacy  and  try  to  make  plausible  and  extend  the  formalism 
used  by  Stern  and  Malkus  (1953)>  by  recourse  to  a physical  argue- 
ment. 

It  is  then  sought  to  test  the  hypothesis  by  applying 
the  theory  to  the  well-known  phenomenon  of  the  sea  breeze.  This 
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will  be  based  on  the  perturbation  of  an  initially  undisturbed 
gradient  wind  that  is  perpendicular  to  a heated  coast.  The  super- 
position  of  this  gradient  wind  is  not  only  a realistic  condition 
that  has  not  been  investigated  theoretically  (to  the  author's 
knowledge),  but  Insures  that  the  linearization  technique  is  sounder 
than  if  it  were  applied  to  an  atmosphere  initially  at  rest.  If 
the  present  theory  of  atmospheric  heating  is  adequate  it  should 
give  predictions  on  the  sea  breeze  component  which  are  in  agree- 
ment with  the  well  known  qualitative  features,  and  in  addition, 
it  may  be  expected  to  lead  to  specific  quantitative  conclusions 
regarding  the  variation  of  the  sea  breeze  with  the  various  para- 
meters, which  are  capable  of  oeing  checked  by  observations. 

A simple  description  of  the  physical  processes  which 
produce  the  mean  motions  may  be  obtained  by  reference  to  Figure  1, 
which  is  a schematic  diagram  of  the  air  flow  over  an  island  that 
is  at  a constant  temperature  above  the  surrounding  water.  At 
large  distances  upstream  the  undisturbed  flow  is  taken  as  adiabatic, 
with  a slight  gradient  wind  (U)  perpendicular  to  the  coast.  The 
undisturbed  lapse  rate  and  stability  are  also  considered  constant. 

As  a result  of  the  large  lapse  rates  that  are  established  over  the 
island,  a turbulent  ground  layer  B is  developed  wherein  heat  is 
transferred  upwards  by  eddies.  This  heating  tends  to  displace  the 
mean  streamlines  upwards;  however,  in  order  to  maintain  the  bound- 
ary condition  of  zero  vertical  velocity  at  the  ground,  there  is  an 
opposing  displacement  due  to  the  stability  of  the  air.  The  nature 
and  behavior  of  this  component  is  quite  similar  to  the  displacements 
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rature  above  the  surrounding  water.  The  islani 
he  y direction  (perpendicular  to  this  section) 
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arising  as  the  result  of  air  flowing  over  a mountain  barrier. 
Whereas  the  turbulent  heating  and  the  displacements  associated 
with  it  are  confined  to  the  ground  layer  B (some  hundreds  of 
meters  in  vertical  extent),  the  mountain  component  is  appreciable 
at  points  far  from  B.  Hence  the  mean  streamlines  outside  the 
ground  layer  are  the  same  as  those  that  would  be  produced  by  an 
equivalent  mountain  whose  shape  depends  on  the  temperature  excess 
of  the  island,  the  undisturbed  wind  and  stability,  and  the  vertical 
extent  of  B. 

II  DERIVATION  OF  THE  PERTURBATION  EQUATIONS 

£ » ^ > £»  tj_  = three  rectangular  coordinates  and  time  (in  c.g.s. 

units)  respectively.  (These  will  be  replaced  later 
by  the  corresponding  dimensionless  coordinates 
x,  y,  z,  t.) 

U,  V = components  of  undisturbed  wind  in  £ , ^ directions, 
u',  v',  w'  = components  of  mean  disturbed  velocities  In^%,'1]  , £ 
directions . 

__  _ ft 

p,  p,  T (or  Tm),  a,  s = 7s — = undisturbed  pressure,  density, 

11  ■‘■m 

absolute  temperature,  lapse  rate,  and  stability, 
p',  p ' , T'  = perturbation  pressure,  density,  and  temperature. 

°1  =TTI  + 0 

f = 2wsinXt  = Coriolis  parameter. 

The  axis  will  be  oriented  so  that  ^ Is  In  the  direction 
of  the  horizontal  temperature  gradient  at  the  ground,  and  the 
4 axis  points  vertically  upwards.  On  performing  a first  order 
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linearized  total  derivative. 
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perturbation  of  the  inviscid  hydrodynamic  equations  and  the 
equation  of  state,  one  obtains  the  following  basic  equations: 


(1)  D-jU*  - fv'  = - | 


(5)  p'  = pRT*  + p'Rf 

Equation  (3)  is  now  used  to  eliminate  the  pressure  from 
(1)  and  (2)  and  equation  (I4.)  is  written  in  terms  of  the  displace- 
ment functions  and<P  . (Subscripts  denote  partial  derivatives.) 

(6)  Dju'^  - fv'^  = jl  pj|  + Dqw^ 

(7)  Dqv'^  + fu'^  = 

(8)  u»  = - U(^  w<  = ; V = U ^ 

Substitute  (8)  in  (6)  and  (7). 

(9)  -DX  + lf^)  - ^ py  + DiTfj, 

(10)  Difj5  - f<V&  + = DlVjl, 

From  (5)  it  follows  that 


It  may  be  shown  (see  Malkus  and  Stern,  1903)  that  the 
second  term  on  the  right  hand  side  leads  to  a small  damping  term 
in  the  final  differential  equation,  which  only  effects  the  ampli- 


tude of  the  displacement  at  large  altitudes.  The  same  is  true 
for  the  terms  involving  , which  have  been  omitted  in  the 

previous  elimination. 

To  investigate  the  influence  of  these  damping  terms 
the  reader  is  referred  to  the  literature  on  the  mountain  wave 
problem  (e.g. , Scorer,  19)4.9 ) . These  terms  are  eliminated  at  the 
outset  of  the  present  discussion  by  writing 


(11)  | TJ 

If  the  flow  were  adiabatic  the  temperature  perturbation 
could  be  eliminated  from  the  above  equations  by  utilizing  the 
fact  that  the  total  derivative  of  the  potential  temperature  (e) 
was  zero.  In  the  present  problem,  however,  an  amount  of  thermal 
energy  is  being  added  at  each  point  which  is  equal  to  the  diver- 
gence of  the  eddy  flux  of  heat.  If  this  quantity  Is  denoted  by 
cp  H » 4»  t-^),  where  cp  is  the  specific  heat  at  constant 
pressure,  then  the  first  law  of  thermodynamics  may  be  written  as 


H = 


1 

9 


de 

dt]_ 


The  first  order  perturbation  leads  to 


(12)  H = DXT»  + w'( P - a)  = DXT»  + U(P  - a)^  . 

Equations  (11)  and  (12)  are  now  used  to  eliminate  p*  and  T»  from  (9). 
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0-3)  DlViV  + + Dl2fi)i;  + " life  7f 

(14)  - - hfa  * + C1  ^4  * 0 


where  T,  the  mean  undisturbed  absolute  temperature  has  been 
written  as  Tm  to  avoid  confusion  with  subsequent  notation,  and 

At  this  point  the  dimensionless  coordinates  (x,  y,  t,  t) 
will  be  introduced,  replacing  t^). 

(1S)  y=r/|l 

f2  Q Q } 

In  addition:  c2  * ; D = +f  “fj)  i 

2 1 2 32  ✓ V 

V = + jT2^  r* 11 


Then  (13)  and  ( lip)  become: 

(16)  dV?  ♦*»  ♦ t>Vy«  + 

(17)  - cjfzz  ~ cl^xy  "^rr  + rfxz  * 0 


With  a constant  value  of  the  Coriolis  parameter  the  elimination 
of  f from  (16)  and  (17)  is  quite  simple.  Since  the  x-axis  was 
chosen  in  the  direction  of  the  horizontal  temperature  gradient 
j/_  may  be  omitted  from  the  final  differential  equation. 

v 


<•&  + Ti^V  V + + 


(16) 
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When  ^ and  H are  written  in  terms  of  Fourier  expansions 
whose  typical  terras  are  e^^x  ei^t'y(e)  and  e*11*  e1^*1  B(z) 
respectively  (henceforth  bars  above  letters  indicate  Fourier 
components),  then  equation  (18)  becomes 

(19)  [o2  - (k+X)2]  j-g  - [k2  - k2  <k+X)2]  y . -^kjig 


The  elimination  of  from  (16)  and  (17)  when  the  lati- 
tudinal change  of  the  Coriolis  parameter  is  considered,  is  more 
difficult.  This  has  been  carried  out  for  the  case  when  the  hori- 
zontal temperature  gradient  along  the  ground  is  in  the  east-west 
direction,  so  that 


> f _ 


= 0 and  = jf  (or  in  dimensionless  units 
= A where  The  result  corresponding  to  (19)  is 


y 

(20) 


c2  - 

1 ' EtAt 


- L“2  ' k2  (kt>)2]  V = .Tl/gt 


III  THE  DESCRIPTION  OF  THE  HEATING  FUNCTION 

If  H and  ^ are  set  equal  to  zero  equations  (19)  and  (20) 
give  solutions  for  the  steady  state  mountain  problem  (Queney,  I9I4.7, 
I9I+8).  However,  the  inhomogeneous  term  (H)  Is  the  essential  driving 
mechanism  which  produces  the  mean  perturbations  over  the  flat 
island,  and  it  must  be  specified  by  another  relation  before  these 
equations  can  be  solved. 

In  order  to  describe  the  turbulent  heating  one  might  be 
tempted  to  use  the  conventional  equation  of  eddy  conduction  and 
set  the  total  derivative  of  the  potential  temperature  equal  to  a 
constant  times  the  second  derivative  of  the  temperature  with 


respect  to  height.  It  is  well  known  that  although  the  eddy 
conductivity  varies  considerably  in  space  and  time,  it  is  often 
possible  to  obtain  a first  orientation  into  many  phenomena  by 
assuming  an  effective  constant  value  for  this  quantity.  It  is 
contended  however,  that  in  the  present  type  of  problem  where  the 
heating  induces  non-negligible  mean  vertical  velocities,  this 
approach  is  inapplicable,  at  least  without  considerable  modifica- 
tion. To  give  a simple  illustration  of  this  inadequacy,  consider 
the  essentially  adiabatic  waves  which  occur  outside  the  ground 
layer,  e.g, , lee  waves.  These  waves  vary  in  magnitude  In  the 
vertical  as  well  as  the  horizontal  direction  and  the  second 
derivative  of  temperature  with  height  is  not  zero.  In  fact, 
application  of  the  simple  eddy  conduction  equation  would  lead  to 
the  untenable  conclusion  that  the  amount  of  heat  added  at  a point 
far  in  the  lee  may  be  comparable  with  the  heat  that  is  supplied 
at  a point  in  the  middle  of  the  turbulent  ground  layer.  If  one  is 
to  use  the  ideas  of  eddy  conduction  at  all,  It  is  necessary  to 
distinguish  between  the  temperature  gradients  that  are  maintained 
by  turbulent  transport  and  those  that  are  due  to  adiabatic  con- 
vective motions.  Tnis  difficulty  might,  in  principle,  be  sur- 
mounted by  using  the  conventional  eddy  equation  with  a conduc- 
tivity that  was  an  order  of  magnitude  less  outside  B than  inside, 
but  because  of  tne  complicated  shape  of  B,  this  formalism  seems 
hopelessly  complex  to  introduce  into  the  hydrodynamic  equations. 

On  the  other  hand,  an  adequate  formalism  must  insure  that  the 
heating  function  approaches  zero  at  far  distances  from  the  island. 
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even  when  mean  vertical  velocities  due  to  the  heating  exist  in 
these  regions  (it  is  assumed,  of  course,  that  there  are  no  other 
mechanisms,  such  as  condensation,  which  are  producing  heat  sources). 
The  following  discussion  is  an  extension  and  elaboration  of  the 
method  of  Stern  and  Malleus  (1953)- 

Because  of  the  fact  that  the  temperature  perturbation 
is  due  to  a combination  of  two  different  mechanisms;  namely, 
turbulent  transport  and  adiabatic  convective  motions,  it  is  in- 
convenient to  attempt  to  describe  the  heating  in  terms  of  this 
quantity  as  the  dependent  variable.  Instead,  the  present  formalism 
evolves  about  the  heating  function  which  is  the  divergence  of  the 
eddy  flux,  and  we  shall  try  to  justify  the  major  premise  that  this 
can  be  determined  independently  of  the  mean  vertical  velocities 
which  it  produces.  This  is  clearly  the  case  at  the  lower  bound- 
ary, as  is  seen  from  Eq.  12,  and  here  the  heating  function  can  be 
determined  from  the  temperature  along  and  the  shape  of  this  bound- 
ary, and  is  not  explicitly  related  to  the  vertical  velocities  that 
are  produced  aloft.  What  can  be  said  about  the  variation  of  the 
heating  function  in  the  vertical? 

In  Figure  1,  consider  the  affect  of  inserting  a series 
of  horizontal  grids,  or  large  screens,  into  the  field  of  motion. 

At  any  point  the  spacing  of  the  grldwork  is  large  compared  with 
the  mean  eddy  size  and  small  compared  to  the  distance  over  which 
one  averages  to  obtain  the  mean  vertical  velocity  and  temperature. 
Then  the  turbulent  eddies  will  pass  through  the  gridwork  relatively 
unaffected  while  the  mean  vertical  velocities  will  be  reduced 
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towards  zero.  The  screens  are  effectively  a solid  barrier  to  the 
mean  motions.  It  is  hypothesized  that  the  eddy  flux  of  heat,  and 
in  particular,  its  divergence,  is  unaffected  by  variations  of  the 
mean  perturbations  due  to  the  imposition  of  these  constraints. 

By  this  means  it  is  possible  to  consider  the  extremely  complex 
turbulent  heating  process  independently;  then,  acting  as  a fixed 
driving  force,  the  heating  function  produces  streamline  displace- 
ments as  determined  by  Equation  18,  and  the  boundary  conditions. 

Accepting  this  hypothesis,  a differential  equation  for 
H is  now  derived  by  applying  the  simple  ideas  of  eddy  conduction 
to  the  model  in  which  the  mean  vertical  velocities  have  been  re- 
duced to  zero.  Denote  the  temperature  in  this  model  by  Te  (£»**]  , 

£,  t)  and  note  that  this  is,  in  general,  different  from  T'  (£,  , 

£,  t),  the  temperature  distribution  in  the  model  whose  mean  motions 
it  is  desired  to  Investigate  (Figure  1).  To  see  this,  consider 
the  region  outside  B where  ftwO.  For  the  original  model  (Figure  1, 
flat  island)  it  follows  from  (12)  that 


w'(P 


0.) 


3 T' 

+ U-jj- 


0 


) 


whereas 


'e 


= 0 and  U- 


+ V- 


w ^r 


TTl 


= 0 


for  the  second  model. 

Applying  the  first  law  of  thermodynamics  to  the  second 
model,  we  obtain  the  following  relations: 


11 


(22)  H = K-^r 


Upon  elimination  of  Te  there  results 


In  deriving  (23)  it  is  assumed  that  K has  a constant  effective 
value  for  the  turbulent  ground  layer  B.  The  fact  that  the  measured 
eddy  conductivity  is  an  order  of  magnitude  less,  far  outside  B is 
of  little  moment  in  this  formalism  because  (23)  insures  that  H is 
small  at  these  points.  It  is  to  be  noted  that  when  the  mean  ver-  • 
tical  velocities  can  be  neglected  (23)  reduces  to  the  familiar 
eddy  conduction  equation,  where  temperature  replaces  H. 

In  addition,  the  boundary  value  of  H follows  directly 
from  the  first  law  and  is  independent  of  any  previous  assumptions 
regarding  the  eddy  transport1. 

(21+)  » °>  tl)  = ( + U^—  + T'(£>7  » °*  tl) 

Returning  to  the  model  in  which  the  horizontal  temperature 
gradient  at  the  ground  is  in  the  direction  of  the  x~axis  and 
writing  (23)  and  (2^)  in  terms  of  the  dimensionless  coordinates 
there  results. 


If  the  ground  is  not  flat  but  has  an  elevation  given  by 
^ > o) , then  one  merely  adds  the  term  U(P  - a)  ^ir 
to  the  right  hand  side  of  (21}.).  ^ 


12  - 


(25)  (^.  + £.)  H = b2 

H(x, o,  t)  = ^gs  (^  + >^)  T'(x,o,t) 

where 

>2  = f^ 

Substituting  the  Fourier  components  for  H and  T*  we  get, 
I(z=o)  = yfe  i(k+A)  T0 

(26)  .*.  IT  = ^gs  i(k+>)  T0  exp  - ^(±  i )*  [k+Al^b’1  zj 
where 


T0  is  the  amplitude  of  the  (k,)0  harmonic  of  the  tempera- 

+xi 

ture  at  the  ground  and  (1  i)  = e— according  as  (k+A)  is 
greater  or  less  than  zero. 


IV  THE  GENERAL  SOLUTION  OF  THE  DIFFERENTIAL  EQUATION  FOR  THE 
MEAN  STREAMLINES 

Upon  substituting  (26)  into  (20)  there  results. 


(27) 


C2  - (k+M21  2 

iff.  [**-**  <*»*]* 

} ' k(ft>)  J * 

Z w 1 

. exp  - £(±  i)^  (k+Al^b-1  z] 
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Solving,  one  obtains 


(28) 

where 


y.-  H(k,^)  £h(z,l^A)  - f(z,k,Jt)J 


(29)  ff(k,X)  = 


k(k+^  [l  - k(k+l)  1 


?o_  k(k+»  jj-grmTJ 

sT®  Qk+A)2  - c2j‘ii^AI  . 


and 


r2  _ k2  -k2_£k+.)>?2.  n _ * ) 

c2  - (k+>)2  ' k(k+A)  } 

f(z,k,)k)  = linear  combination  of  e"rz  and  e+r* 

h(z,k»X ) = exp  - £ (±  i)  Ik+Xl^b*"1  zj 

The  solution  for  the  displacement  Y"  for  an  arbitrary 
temperature  distribution  along  the  ground  is 

(30)  z,t)  = J jl(k,X)  elkx  e1^  £h(z,k,A)  - f(z,k,A^dkdA 

The  solution  is  separated  into  two  parts  and  'Y 2 , where 
yjr  = Yi  - ^2  and 

(31)  tx(x,z,t)  = ^d>ei^t  5T(k,>)  eikx  h(z,k,A)  dk 

/*> 

<*eiAt  / H(k,>)  eikx  f(z,k,A)  dk 


The  first  component  satisfies  the  heat  conduction  equation 

2„ 

' 2 > a il 

(• 


>Vl 

Tt  + 77>  "1  - b “TT5 
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while  the  second  component  satisfies  the  homogeneous  differential 
equation.  Since  ^(x,o,t)  = 0 and  h(o,k,  \)  = 1.  Therefore, 


M(x, t)  has  been  called  the  equivalent  mountain  function 
since  the  component  *»t)  is  mathematically  identical  to  the 

air  flow  over  a mountain  whose  profile  is  given  by  (33).  In 
addition,  (33)  also  determines  the  conduction  component  but 
aside  from  the  fact  that  it  is  necessary  to  satisfy  the  boundary 
condition  at  the  ground  It  is  of  little  interest,  since  It  de- 
creases rapidly  outside  the  ground  layer.  Thus  the  problem  of 
air  flow  over  nested  terrain  is  essentially  reduced  to  the  investi- 
gation of  the  equivalent  mountain  (Equation  33).  In  the  discus- 
sion to  follow  it  Is  convenient  to  consider  the  scale  of  the  heat- 
ing In  a fashion  similar  to  that  used  by  Queney  (I9I48)  in  discuss- 
ing mountain  waves.  For  each  scale  different  sets  of  parameters 
become  Important  and  allow  simplification  of  the  results.  The 
scale  divisions  and  the  approximations  entailed  in  each  are 
summarized  on  the  following  page. 
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Scale 

Size  of  "island" 
L (k  m) 

Simplifications  in 
evaluating  equivalent 
mountain 

Small 

r /jC.  1 

r'/5 

T * 2*  , 1 

> 75 

1 -*1 

Coriolis  parameters 
(c,  ^0)  omitted  from 
equivalent  mountain 

Middle 

c J5  > 71 

k24tl 

omitted 

Large 

* 75 ‘ 

W 

5*  1 

V THE  SMALL  SCALE  SOLUTION  - PREDICTIONS  ON  THE  SEA  BREEZE 

In  this  paper,  only  the  small  scale  solution  will  be 
considered.  The  omission  of  a discussion  of  the  larger  scale 
solution  is  not  primarily  due  to  any  inherent  mathematical 
problems,  butt  rather  to  the  conceptual  difficulties  in  assigning 
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values  for  the  eddy  conductivity  and  the  undisturbed  parameters. 
Moreover,  it  is  felt  that  confidence  in  the  heating  hypothesis 
should  first  be  obtained  by  applying  the  theory  to  such  well 
known  local  phenomena  as  the  sea  breeze  to  test  it3  adequacy. 
However,  formal  manipulation  of  the  larger  3cale  solutions  have 
been  carried  through  and  it  is  believed  that  these  may  be  of 
heuristic  value.  The  author  intends  to  be  able  to  report  on 
this  in  the  near  future. 

For  the  sake  of  definiteness  we  consider  a small  flat 

island  (of  the  order  of  tens  of  kilometers  in  width  and  sensibly 

infinite  in  length)  which  is  at  a uniform  temperature  above  the 

surrounding  water.  The  temporal  variation  of  this  temperature 

difference  is  assumed  to  be  represented  by  a finite  number  of 

Fourier  harmonics,  with  the  predominant  one  being  the  diurnal 

period  (P  = 2lj.  hrs),  then  ^ = ■■  — — (^'-'lO"2).  We  now  make  the 

Pfgs 


small  scale  approximation  and  neglect  ijl  and  Sl  in  comparison  with 
unity.  This  is  valid  for  heat  sources  whose  characteristic 
horizontal  dimension  L is  much  less  than  UP.  Equation  (29)  then 
becomes 


= 


Ik 

& 


1 

+ 1-k2 


If  the  ground  temperature  is  given  by 


r+~ 

TQ(x,  t)  = cos  >vt  / T0(k)  eikx  dk 


then  the  equivalent  mountain  is,  from  (33), 
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<310 


M(x,t)  = 


se  J 


elk*  T„(k) 

15  - k2  + 1 

b2 


dk 


This  differs  from  the  steady  state  solution  obtained  by  Stern 
and  Malkus  (1953)  only  in  the  presence  of  the  cosXt  term.  By 
simple  extension  of  their  results  it  is  easily  shown  that  the 
sea  breeze  component  perpendicular  to  the  windward  shore  of 
the  island  is 


(35) 


where 


. 'T'cosXt  >/5s  cosz 

u = “Tier  *1*2  - »it 


J ~r dp ; 

-<*»* 


x 0 


x = 0 is  at  the  windward  shore 

= t1  1 ^ 

b2 


^cosXt  = temperature  excess  of  the  island 

In  Haurwitz’s  (I9I4.7)  heuristic  model  of  the  sea  breeze  Lt  turns 
out  that  without  friction  the  predicted  sea  breeze  is  90°  out 
of  phase  with  the  temperature.  Defant’s  (1951)  model  of  the 
sea  breeze  shows  a phase  difference  of  I4..7  hours  when  friction 
is  not  considered.  His  theory  while  somewhat  similar  to  the 
present  one  uses  the  conventional  law  of  eddy  conduction  and 
further,  neglects  completely  the  non-linear  inertial  term.  It  is 
Seen  from  the  preceding  analysis  that  the  predicted  sea  breeze 
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component  is  in  phase  with  the  diurnal  temperature  wave,  and  it 
is  not  necessary  to  introduce  friction  to  explain  this  often 
observed  phenomenon.  ^ 

Figure  2 (after  Stern  and  Malkus,  1953)  shows  the  var- 
iation of  the  sea  breeze  component  (u’ ) at  the  windward  shore 
with  the  other  atmospheric  parameters.  The  quantities  u1  and  7" 
in  this  diagram  may  be  interpreted  as  instantaneous  values  in 
the  2l|  hour  cycle.  From  this  diagram  it  is  seen  that  the  ratio 
of  the  eddy  conductivity  to  the  square  of  the  basic  current  is 
of  great  importance  in  determining  the  magnitude  of  the  sea  breeze. 
It  will  be  noted  that,  on  the  scale  considered,  the  Coriolis 
force  has  no  influence  on  u'  and  this  is  due  to  the  fact  that  the 
horizontal  pressure  gradient  and  the  inertial  forces  are  much 
larger  and  control  the  motion  in  the  x-z  plane.  However,  the 
sea  breeze  component  parallel  to  the  coast  is  primarily  balanced 
by  the  Coriolis  force,  as  shall  be  seen  subsequently  when  the  hod- 
ograph  Is  discussed.  Before  turning  to  this,  some  observational 
work  in  connection  with  the  preceding  formula  will  be  discussed. 

The  success  of  the  present  model  in  obtaining  a 
realistic  picture  of  the  sea  breeze  phenomenon  has  encouraged 
attempts  to  obtain  quantitative  observational  checks.  By 
differentiating  equation  (35)  with  respect  to  x and  denoting 
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simple  expression: 


where 

= temperature  excess  of  island  at  a given  time, 

U = undisturbed  wind  component  perpendicular  to  coast, 

Tm  = mean  absolute  atmospheric  temperature. 

Although  all  explicit  reference  to  the  vertical  distribution  of 

heating  has  been  eliminated  in  the  above  equation,  its  validity 

depends  on  the  fact  that  a constant  effective  eddy  conductivity 

for  B may  be  chosen.  In  addition  it  contains  relatively  simple 
✓ 

observables  (temperature  and  horizontal  wind  speed)  and  seems 
capable  of  being  checked.  Its  validity  would  not  only  confirm 
the  theory,  but  would  afford  a simple  means  of  calculating  K 
(see  Stern  and  Malleus,  1953) » by  means  of  an  additional  formula. 

During  the  summer  of  1952,  a small  exploratory  field 
program  was  undertaken,  mainly  with  the  purpose  of  determining 
whether  or  not  the  expected  slope  in  the  wind  profile  indicated 
by  (36)  could  be  measured.  Figure  3 shows  a series  of  wind 
profiles  obtained  running  from  the  southern  shore  of  Nantucket 
to  the  interior  of  the  Island.  Three  minute  averages  on  a small 
cup  anemometer  placed  about  seven  feet  above  the  ground  were 
made  at  each  station.  Although  measurements  of  distance  and 
angles  were  rough,  determination  of  the  mean  slope  for  each  of 
the  runs  was  made  and  the  quantity  computed.  The  results 

are  shown  in  Table  1.  It  is  to  be  pointed  out  that  for  comparison 


^ Observed,  horizontal  profiles  of  the  sea  breeze  at  the  windward 

of  Nantucket  Island. 
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with  the  theoretical  relation  (36),  what  is  really  wanted  Is  the 
wind  gradient  from  the  shore  out  to  sea  rather  than  the  one  that 
was  measured.  Table  1 indicates  that  the  measured  quantity  ^ 
is  constant  within  observational  accuracy,  fluctuating  about  a 
mean  of  two.  The  systematic  deviation  from  the  predicted  value 
of  unity  is  believed  to  be  due  to: 

a)  The  ratio  of  the  measured  slope  (on  land  where  the  sea 
breeze  is  decreasing  with  x)  to  the  slope  over  the  sea 
(where  the  sea  breeze  is  increasing). 

b)  The  non-uniformity  of  the  island  temperature. 

A more  intensive  observational  program  is  being  considered  by 
the  meteorology  group  at  this  Institution  which  is  based  on  the 
measurement  of  pressure  differences  to  test  (36).  Thus,  since 


' equation  (36)  becomes 

V P V 

= 1 


where 

p = mean  atmospheric  density  1.2  x 10"3  gms  cm"3 
g = 980  cm  sec”2 

= temperature  excess  of  island  at  a given  time 

O _ | 

- -y~~-  = W = perturbed  pressure  gradient  at  windward  shore  at  a 
given  time. 

There  seems  to  be  little  doubt  that  the  pressure  gradient  associated 
with  the  sea  breeze  can  be  accurately  measured  (see  for  example, 
Leopold  (I9I4.9)). 


Table  I 


RESULTS  OF  SEA  BREEZE  MEASUREMENTS 


Run  No. 

G 

cm/sec/800  ft. 

Y 

°F 

Angular  Deviation 
of  undisturbed 
wind  from  normal 
to  coast 

1 

95 

3.7 

45° 

1.8 

2 

90 

5.0 

0 

0 

ro 

2.2 

3 

100 

5.0 

20° 

1+ 

60 

5.7 

0° 

Computation  of  YiL  from  observational  measurements  for  com- 
■rtGTjnO 

parison  with  theoretical  results.  U Is  the  component  of  un- 
disturbed gradient  wind  that  is  perpendicular  to  the  coast, 

'T"  is  the  temperature  excess  of  the  island  above  the  water,  G 
is  the  horizontal  gradient  of  the  horizontal  wind,  Tm  = 300°C, 
g = gravity,  it  = 3.1I4.I, 
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In  concluding  this  section  we  present  some  conclusions 
on  the  turning  and  the  hodograph  of  the  sea  breeze.  This  subject 
has  received  considerable  theoretical  and  observational  attention 
and  it  is  generally  accepted  that  the  turning  is  due  to  the 
earth's  rotation.  The  major  purpose  of  the  following  development 
is  to  show  how  the  various  atmospheric  parameters  such  as 
stability,  eddy  conductivity  and  undisturbed  wind  speed  effect 
the  hodograph.  To  compute  v'  equation  (2)  is  used^ and  since  the 
temperature  gradient  in  the  y direction  vanishes,  the  y derivatives 
of  the  perturbation  quantities  are  neglected.  In  dimensionless 
coordinates  the  y equation  of  motion  becomes 


<Tt  + Tx'  v'  = " c 


The  solution  of  this  equation  is 

v'  (x,t)  = - c jT  u'  £x  - t +J  d 7J 

where  a is  a constant  of  integration,  as  may  be  verified  by 
differentiation.  This  will  be  evaluated  at  x = 0,  which  is 
the  windward  shore  of  the  island.  Now,  u'  varies  as  cos^t, 
hence  at  t = - ^ , u'  = 0,  and  we  assume  v'  = 0 at  t = tQ  - 
where  tQ  is  small  compared  to  ^ . 

(37)  ..  „ /u.  <7,  J ♦ t)  dj 

tQ-  £ -t 


The  value  of  u'  is  given  by  equation  (3£),  but  rather  than 
attempt  to  integrate  this  complicated  expression  it  is  approximated 
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by  a simple  exponential  function  which  fits  it  in  the  vicinity  of 
the  windward  shore.  Thus,  if  u0  is  the  amplitude  of  u'  given  by 
(35)  when  x * t = 0 and  G is  the  slope  of  this  curve  at  the  same 

point  ( i. e. , G * If  ) then, 

u’(x,  t)»u0  exp  cos>t  x<0 

*0\ 


where 


u°  * £ ! $ ln  KTgr 
G = 


We  seek  the  behavior  of  v'  when  t>^  - h”  •*  10  p i • 0 • 9 well 

after  the  time  of  onset  of  the  sea  breeze.  Since  x A 0 

UOi 


Ui 


when  x = tQ  - -?£ — t we  may  write  (3?)  as 

41  ♦ u7, 


^ * /* 


e*^  ©xp  (- 


(38)  ’ * 

‘ '*  -C  U0 


Re 


,3>t 


+ a 


Uo 


where  Re  denotes  "real  part  of". 

Prom  equation  (38)  it  is  seen  that  the  sea  breeze  component 
parallel  to  the  coast-  is,  in  general,  not  of  the  same  phase 


as  the  diurnal  temperature  wave.  The  quantity 


“op 


is  a 


measure  of  the  deoth  of  the  sea  breeze,  or  how  far  out  to  sea 
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it  extends.  When  the  depth  is  small  compared  with  the  distance 
that  an  air  parcel  travelling  with  the  speed  of  the  undisturbed 
wind  would  cover  in  hours,  then  v*  = -c  u0  cosXt.  In  this  case 
v'  is  much  smaller  (~10" 2 times)  than  u*  and  the  hodograph  of  the 
sea  breeze  components  would  be  a straight  line  almost  perpendicu- 
lar to  the  coast.  Figure  1+  shows  the  theoretical  hodograph 
when Q = .05*  X=  c = 10"2.  This  might  correspond  to  a 

sea  breeze  that  was  five  kilometers  deep.  It  Is  readily  shown 
from  (38)  that  these  hodographs  are  elliptic.  The  point  labelled 
At  = 0 In  Figure  I4.  represents  the  time  of  maximum  heating,  I.e., 
the  temperature  difference  between  land  and  sea  Is  a maximum. 

For  deeper  sea  breezes,  - decreases  thereby  increasing  both 


Uq 


7u  2 


the  amplitude  and  phase  angle  of  v’,  and  hence  the  hodograph 
would  tend  to  become  more  circular.  These  relations  between  the 
depth  of  the  sea  breeze  and  the  shape  of  the  hodograph  at  the 
coast  would  appear  to  have  some  interesting  practical  applica- 
ti ons . 


SUMMARY  AND  CONCLUSIONS 

The  first  part  of  this  paper  deduces  the  equation  of 
motion  for  a non-adiabatic  atmosphere  where  the  mean  motions 
are  small  perturbations  compared  with  the  velocities  in  the  basic 
current.  The  resulting  partial  differential  equation  is,  aside 
from  the  non-homogeneous  forcing  function,  the  same  as  that  for 
air  flowing  over  a mountain. 
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It  is  then  sought  to  determine  the  mean  perturbations 
that  have  their  energy  source  in  the  turbulent  heating  near  the 
ground.  The  inadequacy  of  the  conventional  formalism  which  de- 
scribes the  heating  by  means  of  an  eddy  conduction  equation  for 
the  potential  temperature  is  shown.  Instead  the  heating  is  de- 
scribed in  terms  of  H which  is  equal  to  the  divergence  of  the 
eddy  flux  of  heat,  or  to  first  order  the  total  derivative  of  the 
potential  temperature.  On  the  basis  of  some  physical  arguements 
and  hypotheses  it  is  proposed  that  this  function  satisfies  the 
equation 


where  _SL  is  the  total  hydrodynamic  derivative  and  K is  a mean 
dt 

value  of  the  eddy  conductivity  for  the  turbulent  ground  layer. 

The  steady  state  solution  of  the  equations  of  motion 
using  a heating  function  that  satisfies  the  above  equations  has 
been  discussed  by  Stern  and  Malkus  (19$3)  and  compared  with 
observations  over  Nantucket  Island. 

In  order  to  further  substantiate  the  theory  this  paper 
then  proceeds  to  discuss  the  sea  breeze  by  retaining  the  time 
derivatives  and  the  Coriolis  parameter  where  necessary  on  the 
scale  pertaining  to  this  phenomenon.  It  is  shown  that  the  sea 
breeze  component  perpendicular  to  the  coast  is  in  phase  with 
the  diurnal  temperature  wave  and  that  it  is  not  necessary  to 
Introduce  friction  to  explain  this.  However,  the  sea  breeze 
component  parallel  to  the  coast  has  a variable  phase  depending 


1 
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upon  the  depth  of  the  sea  breeze.  The  theory  demonstrates  the 
role  of  the  various  atmospheric  parameters  on  the  shape  and  size 
of  the  elliptic  hodograph  of  the  sea  breeze  components. 

Preliminary  observations  of  the  sea  breeze  produced 
by  a small  flat  island  (Nantucket,  Massachusetts)  have  indicated 
that  quantitative  as  well  as  qualitative  agreement  of  the  theory 
may  be  expected.  By  means  of  this  theory  the  sea  breeze  may  be 
used  as  a fruitful  tool  to  investigate  the  gross  properties  of 
turbulent  heating  near  the  ground,  and  it  is  hoped  to  be  able 
-to  continue  these  observational  studies. 

The  satisfactory  picture  of  local  convective  phenomena 
that  is  produced  by  turbulent  heating  near  the  ground,  would 
suggest  that  the  theory  now  be  applied  to  a larger  scale.  Pre- 
liminary theoretical  investigations  indicate  that,  despite 
obvious  difficulties  arising  from  assigning  eddy  conductivities 
on  this  scale,  certain  conclusions  can  be  drawn  which  may  be  of 
heuristic  value.  This  will  constitute  the  subject  of  a future 
paper. 
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Spring  19,  Maryland  1 
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Bureau  of  Ships , Navy  Department,  Washington  2$,  D.  C. 

Attention:  Code  8$1  1 

Bureau  of  Ships,  Navy  Department,  Washington  2$,  D.  C. 

Attention:  Code  8II4.  1 

Bureau  of  Ships,  Navy  Department,  Washington  2$,  D.  C. 

Attention:  Code  327  2 

Chief  of  Naval  Operations,  Navy  Department,  Washington  2$, 

D.  C.  Attention:  OP-533^  2 

Oceanographic  Division,  0.  S.  Navy  Hydrographic  Office, 

Suit land,  Maryland  1 

Library,  Naval  Ordnance  Test  Station,  Inyokern,  China 

Lake,  California  1 

Project  AROWA,  U.  S.  Naval  Air  Station,  Building  R-I4.8, 

Norfolk,  Virginia  1 

The  Chief,  Armed  Forces  Special  Weapons  Project,  P.  0. 

Box  2610,  Washington,  D.  C.  1 

Office  of  the  Chief  Signal  Officer,  Engineering  and 
Technical  Service,  Washington  2£,  D,  C.  Attention: 

SIGGE-M  1 

Meteorological  Branch,  Evans  Signal  Laboratory*  Belmar, 

New  Jersey  1 

Office  of  the  Quartermaster  General,  2nd  and  T Streets, 

S.  W. , Washington  2$,  D.  C.  Attention:  Environmental 

Protection  Section  1 

Office  of  the  Chief,  Chemical  Corps,  Research  and 
Engineering  Division,  Research  Branch,  Army  Chemical 
Center,  Maryland  2 

Commanding  Officer,  Air  Force  Cambridge  Research  Labora- 
tories, 230  Albany  Street,  Cambridge,  Massachusetts 
Attention:  ERHS-1  1 

Chief  of  Staff,  Headquarters  USAF,  The  Pentagon,  Washington 
2$,  D.  C.  Attention:  AFDRD-RE  1 
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Headquarters * Air  Weather  Service , Andrew  A.  F.  Base* 

Washington  20*  D.  C.  Attention:  Director  Scientific 

Services  2 

Commanding  General*  -Air  Materiel  Command*  Wright- Patterson 
Air  Force  Base*  Dayton*  Ohio  Attention:  MCREEO  1 

Commanding  General*  Air  Force  Cambridge  Research  Center* 

230  Albany  Street*  Cambridge  39*  Massachusetts 

Attention:  CRHSL  1 

Commanding  General*  Air  Research  and  Development  Command* 

F.  0.  Box  1395*  Baltimore  3*  Maryland  Attention:  RDDG  1 

Department  of  Meteorology*  Massachusetts  Institute  of 
Technology*  Cambridge*  Massachusetts  Attention:  H. 

G.  Houghton  1 

Department  of  Meteorology*  University  of  Chicago*  Chicago 

37*  Illinois  Attention:  H.  R.  Byers  2 

Institute  for  Advanced  Study*  Princeton,  New  Jersey 

Attention:  J.  von  Neuman  1 

Scripps  Institution  of  Oceanography,  La  Jolla*  California 

Attention:  R.  Revelle  1 

General  Electric  Research  Laboratory*  Schenectady*  New 

York  Attention:  I.  Langmuir  1 

St.  Louis  University,  3621  Olive  Street,  St.  Louis  8* 

Missouri  Attention:  J.  B.  Macelwane*  S.  J.  1 

Department  of  Meteorology,  University  of  California  at 
Los  Angeles*  Los  Angeles*  California  Attention: 

M.  Neiburger  1 

Department  of  Engineering,  University  of  California  at 
Los  Angeles,  Los  Angeles*  California  Attention: 

L.  M.  K.  Boelter  1 

Department  of  Meteorology*  Florida  State  University* 

Tallahassee*  Florida  Attention:  W.  A.  Baum  1 
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Woods  Hole  Oceanographic  Institution!  Woods  Hole, 

Massachusetts  Attention:  C.  Iselin  1 

The  Johns  Hopkins  University,  Department  of  Civil 
Engineering.  Baltimore,  Maryland  Attention:  R.  Long  1 

The  Lamont  Geological  Observatory,  Torrey  Cliff, 

Palisades,  New  York  Attention:  M.  Ewing  1 

The  Johns  Hopkins  University,  Department  of  Physics, 

Homewood  Campus,  Baltimore,  Maryland  Attention: 

G.  Plass  1 

New  Mexico  Institute  of  Mining  and  Technology,  Research 
and  Development  Division,  Socorro,  New  Mexico 
Attention:  E.  Workman  1 

University  of  Chicago,  Department  of  Meteorology, 

Chicago  37,  Illinois  Attention:  H.  Riehl  1 

Woods  Hole  Oceanographic  Institution,  Woods  Hole, 

Massachusetts  Attention:  A.  Woodcock  1 

Iowa  State  College,  Department  of  Physics,  Ames,  Iowa 

Attention:  J.  E.  MacDonald  1 

Blue  Hill  Meteorological  Observatory,  Harvard  University, 
Milton  86,  Massachusetts  Attention:  C.  Brooks  1 

Department  of  Meteorology,  University  of  Washington, 

Seattle  $,  Washington  Attention:  P.  E.  Church  1 

Laboratory  of  Climatology,  Johns  Hopkins  University, 

Seabrook,  New  Jersey,.  Attention:  C.  W.  Thornthwaite  1 

Institute  of  Geophysics,  University  of  California  at 
Los  Angeles,  Los  Angeles,  California  Attention: 

J.  Kaplan  1 

Department  of  Meteorology,  New  York  University,  New 

York  53*  New  York  Attention:  B.  Haurwits  1 

Texas  A and  M,  Department  of  Oceanography,  College 

Station,  Texas  Attention:  D.  Leipper  1 
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Massachusetts  Institute  of  Technology,  Department  of 
Meteorology,  77  Massachusetts  Avenue,  Cambridge, 

Massachusetts  Attention:  T.  F.  Malone  1 

Cornell  University,  Department  of  Agronomy,  Division 

of  Meteorology,  Ithaca,  New  York  1 

Pennsylvania  State  College,  School  of  Mineral  Indus- 
tries, Department  of  Earth  Science,  State  College, 

Pennsylvania  Attention:  H.  Neuberger  1 

Rutgers  University,  College  of  Agriculture,  Department 

of  Meteorology,  New  Brunswick,  New  Jersey  1 

University  of  Texas,  Department  of  Aeronautical 

Engineering,  Austin,  Texas  Attention:  K.  H.  John  1 

University  of  Utah,  Department  of  Meteorology,  Salt 

Lake  City,  Utah  Attention:  V.  Hales  1 

University  of  Wisconsin,  Department  of  Meteorology, 

Madison,  Wisconsin.  Attention:  V.  Suomi  1 

National  Advisory  Committee  of  Aeronautics,  172l|.  F Street, 

N.  W. , Washington  2$.  D.  C.  2 

U.  S.  Weather  Bureau,  2ij.th  and  M Street,  N.  W. , 

Washington  25,  D.  C.  Attention:  Scientific  Services 

Division  2 

Committee  on  Geophysics  and  Geography,  Research  and 

Development  Board,  Washington  25,  D.  C.  2 

Air  Coordindating  Committee,  Subcommittee  on  Aviation.. 
Meteorology,  Room  2D889-A,  The  Pentagon,  Washington  25, 

D.  C.  1 

American  Meteorological  Society,  3 Joy  Street,  Boston 

8,  Massachusetts,  Attention:  The  Executive  Secretary  1 

Laboratory  of  Climatology,  Johns  Hopkins  University, 

Seabrook,  New  Jersey  Attention:  M.  Halstead  1 

U.  S.  Naval  War  College,  Newport,  Rhode  Island  1 
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Brookhaven  National  Laboratory,  Upton,  L.  I.,  N.  Y. 

Attention:  Meteorology  Group  1 

Chief,  Meteorological  Division,  Biological  Department 

Chemical  Corps,  Camp  De trick,  Frederick,  Maryland  1 

Dr.  August  Raspet,  Engineering  and  Industrial  Research 
Station,  Mississippi  State  College,  State  College, 
Mississippi  2 

Dr.  E.  W.  Hewson,  Diffusion  Project,  Round  Hill,  South 

Dartmouth,  Massachusetts  1 

Dr.  Hunter  Rouse,  Director,  Iowa  Institute  of  Hydraulic 

Research,  State  University  of  Iowa,  Iowa  City,  Iowa  1 

Head,  Department  of  Physics,  University  of  New  Mexico, 

Albuquerque,  New  Mexico  1 

Mr.  Wendell  A.  Mordy,  Hawaiian  Pineapple  Research 

Institute,  Honolulu,  Hawaii  1 


